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ABSTRACT 

We present the results of a search for high-energy 7-ray emission from a large sample 
of galaxy clusters sharing the properties of three existing Fermi-hAT detections (in 
Perseus, Virgo and Abell3392), namely a powerful radio source within their brightest 
cluster galaxy (BCG). From a parent, X-ray flux- limited sample of 863 clusters, we 
select 114 systems with a core-dominated BCG radio flux above 50/75 mJy (in the 
NVSS and SUMSS surveys, respectively), stacking data from the first 45 months of 
. the Fermi mission in three energy bands, to determine statistical limits on the 7-ray 

O f fluxes of the ensemble of candidate sources. 

For a >300 MeV selection, the distribution of detection significance across the 
sample is consistent with that across control samples for significances <3ct, but has a 
tail extending to higher values, including three >4cr signals which are not associated 
with previously identified 7-ray emission. Modelling of the data in these fields results 
in the detection of four non-2FGL Fermi sources, though none of these appear to be 
unambiguously associated with the BCG candidate. Only one is sufficiently close to 
be a plausible counterpart (RXC J0132.6— 0804), and the remaining three appear to 
be background AGN. A search at energies >3 GeV hints at emission from the BCG 
in A 2055, which hosts a BL Lac object. 

There is no evidence for a signal in the stacked data, and the upper limit de- 
rived on the 7-ray flux of an average radio-bright BCG in each band is at least 
an order-of-magnitude more constraining than that calculated for individual objects. 
-PiGev/-Fi.4GHz for an average BCG in the sample is <15, compared with «120 for 
NGC 1275 in Perseus, which might indicate a special case for those objects detected at 
high energies. The tentative suggestion that point-like beamed emission from member 
galaxies comprise the dominant bright 7-ray sources in clusters, implies searches for 
evidence of dark matter annihilation or large scale merger shock signatures for exam- 
ple, need to account for a significant level of contamination from within each cluster, 
that is both highly stochastic and varies significantly over time. 

Key words: galaxies: active - galaxies: clusters: general - gamma rays: observations 
- radio continuum: general - radiation mechanisms: non-thermal - (galaxies:) cooling 
flows 



1 INTRODUCTION 

Feedback from the central active galactic nucleus (hereafter 
AGN) within a cooling-core cluster of galaxies is consid- 
ered to play an integral large-scale role in the suppression 
of the observed cooling flow (see e .g. I Fabian et all |2002| ; 
ISiiacki et af]|2007l : lBower et al-lhoOol h The hot, diffuse, intr- 
acluster gas, which comprises 90% of the baryonic matter of 
the cluster and is bound gravitationally by a cold dark mat- 
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ter (DM) halo, is revealed observationally through thermal 
Bremsstrahlung emission in the X-ray band. From velocity 
dispersion measurements, the atmosphere of this intraclus- 
ter mediu m (ICM) typically has a tem perature between 10 7 
and 10 s K (|McNamara fc Nulsenl2007l 'l. The gas loses energy 
to radiation primarily in the X-ray band, condensing onto 
the central galaxy, and sub-sonic, pressure-driven inflows de- 
velop within the dense cooling region, where the radiative 
cooling tim e for such em ission t coo \, is less than the age of 
the system (|Fabianlll99"i 'l. 

However, as evidenced by spectroscopic data from 
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current-generation X-ray observatories Chandra and XMM- 
Newton, emission predicted by the standard cooling-flow 
model from expected quant ities of cooling gas is not de- 
tected (|Peterson et al.ll2003h . This soft X-ray deficit neces- 
sitates some reheating mechanism that acts to replenish the 
energy being radiated away, thus quenching the cooling flow. 
An attractive agent for this feedback is the central AGN, 
whose outbursts would not only resolve the dearth of X- 
ray luminosity, but also account for the measured trunca- 
tion at the high-ma ss end of the galaxy luminosity function 
jBenson et alj|2003l ). 

Support for this explanation is given by high-resolution 
obse rvations of cavities in the X-ray halo of m any systems 
(e.g. iFabian et al. l l200d : lMcNamara et alJboOCh . coincident 
with radio synchrotron emission associated with the activ- 
ity of the central engine. The AGN in the brightest clus- 
ter galaxy (hereafter BCG) injects energy into the ICM in 
the form of bubbles containing relativistic particles, which 
displace the thermal gas. So-called 'ghost cavities' have also 
been found using Chandra, interpreted as relics of past AGN 
activity; bubbles that have previously detached from the ra- 
dio lobes and risen buoyantly through the ICM, no longer 
bright at 1.4 GHz. They are now known to be fille d with 
low-frequency radio sy nchrotron emission (e.g. IClarke et al.l 
120051 ; I Wise et alj|2007h . This evidence implicates AGN out- 
bursts as intermittent phenomena; the power output aver- 
aged over time of the central engine able to balance the ra- 
diative cooling of the cluster c ore. For a review of AGN heat - 
ing in clusters of galaxies, see iMcNamara fc Nulsenl (|2007l ). 

Both models of galaxy formation/clustering, and of 
feedback, and observations (including those aforementioned 
in the hard X-ray and radio bands), establish galaxy clus- 
ters a s hosts to significan t non-thermal pa rticle popula- 
tions (|Volk fc Atovanll 19991 ; felasi et al.ll200Sf ). The remain- 
ing tracer of cosmic ray (CR) acceleration is 7 radiation, 
and several emission scenarios predict BCGs to be sources 
of high-energy (HE) and very-high-energy (VHE) photons, 
with spectra dependent on the dominant energy content of 
the AGN jets; the nature of the seed particles they inject into 
the ICM, the strength of the local magnetic field, and the 
targe t density around the inflated bubbles fe.g. lHinton et ail 
12007ft . 

Ultra-relativistic particles may play an important role 
on several different scales in and around BCGs in cooling- 
core clusters, for example in AGN jets, mini-halos, and 
cluster-scale halos. 7-ray emission m ay arise in each via 
a nu mber of dissimilar processes (e.g. IPfrommer fc E nfilin 
120041 ). Primary CR electrons accelerated in jets, or re- 
accelerated fossil CRs in mini-halos will produce inverse 
Compton (IC) emission. Inelastic proton-proton collisions 
of hadronic CRs on the target medium will lead to 7r de- 
cay e mission and IC emission from secondary (or ca scade) 
(e.g. iBlasi fc Colafrancescol 1 19991 ; lAharonianl 120021 ) elec- 
trons. Most of these BCG-driven processes will produce GeV 
emission that is point-like for current detectors. However, it 
should be noted that - being massive and DM-dominated 
- clusters of galaxies are expect ed sources of extended 7 - 
ray emission via DM annihilation (|Ackermann et aI1l2010al ). 
Similarly, treating galaxy clusters as long-term reservoirs of 
CR hadrons, diffuse HE emission from ir° decay of particles 
accelerated (f or example) in cluster-scale accretion shocks is 
theorised fe.g. lPfrommer fc Enfflinll2004l ). The various emis- 



sion scenarios produce different spectral, temporal and spa- 
tial signatures at GeV energies. The repository of HE data 
afforded by the Fermi Gamma-ray Space Telescope (here- 
after Fermi) since its launch in 2008, is therefore an at- 
tractive resource when used in conjunction with multiwave- 
length data. Whilst a detection of extended cluster emissio n 
at GeV energies is yet to be made l|Ackermann et alll2010bl ) , 
point-like emission from a number of BCGs is seen and there 
is considerable potential for improving our understanding of 
this class of objects using GeV data. 

To date, excluding systems whose AGN jets are 
favourably aligned close to the observer's line of sight so 
as to boost emission via relativistic beaming - associated 
with 'BL Lac' objects w ithin the unified model of active 
galaxies (| Antonuccill 19931 ) - only two BCGs have been con- 
firmed as emitters of HE 7 rays. These are NGC 1275 and 
M87 (|Abdo et al.ll2009al lbl): the dominant member galaxies 
of the cooling-core clusters Perseus and Virgo, respectively. 
Both are known to host bright ce ntral radio sources asso - 
ciated with the act ive n ucleus (see Vermeulen et al. Il994l); 
lAsada et al. I (|2006t ) and iBiretta et al. I jl99lMSparks et al.l 
Jl996t)1. and are also detected at ve ry-high energies 
|Aleksic et al]|2012l ; lAharoman et alj|2006j ). The multiwave- 
length emission of AGN is characteristically variable, and 
indeed the 7-ray flux is found to vary on timescales dow n 
to days in the cas e of NGC 1275 (uBrown fc Adams! |201 il l; 
similarly for M 87 (|Abramowski et alj|2012l ). Light-crossing- 
time arguments thus restrict the physical size of the emitting 
region responsible for, at least a component of, the 7-ray 
emission. 

These sources are HE laboratories for the class of BCGs, 
and motivate a search for further instances of 7-ray activ- 
ity within it. The purpose of this work is to conduct such 
a search, using data from the Large Area Telescope (LAT) 
on board the Fermi satellite, accumulated thus far over the 
mission lifetime. To this end, we set out to imitate a deeper 
observation than is currently achievable in HE 7 rays given 
the sensitivity of the Fermi-LAT, by constructing a sample 
of many potential sources within the BCG class, and stack- 
ing the output of individual analyses. 

We defend our derivation of a suitable sample of candi- 
date sources in [J2] and the Fermi-LAT analysis procedure 
then carried out for each target is described in SJ31 The re- 
sults are presented in Sj4] and discussed in Sj5] Conclusions 
of the work are then drawn in Sj6] 



2 TARGET SELECTION 

Considering NGC 1275 to be prototypical of 7-ray sources 
within the class of BCGs, our selection criteria were there- 
fore designed to identify similar objects. The key require- 
ments were the presence of a massive cluster of galaxies 
and evidence of non-thermal activity in the BCG. We con- 
structed a sample of candidates drawn from five parent cata- 
logues of X-ra y flux-limited gala xy clusters: the Brightest 55 
(B55) Sample (lEdge et al.lll99Cl). the ROSAT Bright Cluster 
Sample fBCS) labeling et al.lll99sT), the extended RO SAT 
Bright Cluster Sample feBCS) labeling et alj 120001 ). the 
ROSAT-ESO Flux-Li mited X-ray (REFLEX ) Galaxy Clus- 
ter Survey Catalogue l|Bohringer et al. 20041) , and the MAs- 
sive Cluster Survey (MACS) labeling et all l200lf ). encom- 
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passing objects from both the Southern and Northern hemi- 
spheres and including the brightest clusters from low Galac- 
tic latitudes. We include two clusters (RXC J1350.3+0940 
and RXC J1832. 5+6848) that meet the eBCS flux limit 
but were mis-identified as AGN when the sample was pub- 
lished, due to the presence of a bright radio source, and 
for comparison, two X-ra y-luminous clusters that contai n 
known AGN (4C +55.16 dHlavacek-Larrondo et al.l 1201 J) , 
and E 1821+643 l|Russell et all l20ld )1 though the cluster 
emission falls just below the eBCS limit once the AGN con- 
tribution is taken into account. 

From this sample of 863 unique clusters (accounting for 
the small overlap between these five samples), we selected 
the 151 systems in which the BCG is de tected above 50 mjy 
in the NRAO VLA Sky Survey (NVSS) i|Condon et al.lll998T ) 
or 75 mjy in the Sydney Un iversity Molonglo Sky Survey 
(SUMSS) l|Mauch et al.H2003T ). We took care to remove sys- 
tems where the radio emission is from a cluster member 
projected close to the BCG. Of this sample we have higher 
resolution radio imaging from the VLA or ATCA telescopes 
at 5 and/or 8 GHz for the majority (114 of 151), which can 
be used to identify tho se sources with a s ignificant, flat spec- 
trum core component (|Hogan et alJlprepT) . 

We have optical spectra for 148 of the 151 BCGs recog- 
nised by our radio selection, and of these 64% (97 objects), 
exhibit optical line emission. This is significantly higher 
than the 25-30% of BCGs in an X-ray selected sample 
presenting line emission, indicating that our radio selec- 
tion strongly prefers clusters with a cooling flow (|Burn£ 
19901 ; ICavagnolo et al.ll2008l ; iMittal et al1l2009l : iHogan et all 
pred) . 

The sample represents a complete selection of radio- 
bright clusters that is dominated by cooling-core clusters, 
spanning a wide range in redshift (0.009 < z < 0.45). It can 
be divided into BCGs whose radio emission has some con- 
tribution from the core or is core-dominated, and those with 
extended radio morphology and a lobe-dominated radio flux. 
The former comprise 114 objects and constitute the 'main 
sample' (see Figure Q] and Table [TJ, with the remaining 37 
BCGs forming an 'extended sample'. We compiled a counter 
sample of 65 'control' clusters drawn from the parent X-ray 
sample that do not contain a radio-bright BCG, nor optical 
emission lines, but share the same redshift distribution as 
the radio bright sample for a representative comparison. We 
also generated a pair of control samples of 114 candidate 
sources each with random positions on the sky (generated 
by performing first a positive and then a negative 5° shift 
in Galactic latitude on the coordinates of the main sample 
candidates) for direct comparison with the main sample. 

The parent clusters of the main sample BCGs include 
hosts to five BL Lac objects, two of which are known HE j- 
ray sour ces listed in the Fe rmi-LAT Second Source Catalog 
(2FGL) |Nolan et alj|2012j ): 2FGL J0627. 1—3528 in A 3392 
and 2FGLJ1958.4-3012 in RXC J1958.2— 3011. The other 
three BL Lacs, RGB J1144+676 in A1366, RGBJ1518+062 
in A2055 and B2 2334+23 in A 2627, are not known Fermi 
sources. These BL Lacs are included as the original clus- 
ter identification from the ROSAT Survey incorrectly at- 
tributed the (bulk o f the) X-ray emission t o the cluster (e.g. 
RXC J1958. 2-3011. [BShringer et al1l2004T l. However, in all 
five cases, pointed Chandra or XMM-Newton observations 
indicate that the BL Lac dominates the emission, so the 



cluster does not meet the X-ray flux limit applied to se- 
lect the other clusters. Therefore even this small Fermi-LAT 
detection fraction is an upper limit, given only NGC 1275 
(a misaligned blazar) and M87 (a radio galaxy) are domi- 
nated by the extended X-ray emission from their host clus- 
ter, and hence meet the cluster selection criteria. For the 
purpose of the analyses detailed hereafter, these two well- 
studied BCGs are not included in the sample, since our aim 
is to identify whether any additional sources in the class 
are detected using Fermi . There is no clear cluster around 
RXC J1958.2— 3011, so its classification as a BCG is not cer- 
tain, but with this caveat in mind it is kept, particularly as 
an additional Fermi source to check the output of analyses 
against. 



3 Fermi-LAT OBSERVATIONS AND DATA 
ANALYSIS 

The Fermi-LAT is a pair conversion instrument sensitive to 
photons of energy ~20 MeV - >300 GeV, with a wide field of 
view (2.4 sr) and large effective area («8000 cm 2 on axis at 
1 GeV ). The detector is described in detail in lAtwood et al.l 
(2009) and references therein. Fermi observations are taken 
primarily in all-sky-survey mode, in which the spacecraft is 
pointed at some angle from the zenith, and rocked north and 
south of its orbital plane, scanning the whole sky every ~3 
hours. 

The 151 BCGs of our main and extended samples were 
treated as candidate point sources of HE 7 rays, and each 
analysed using the Fermi Science Tools version v9r27pl, 
and instrument response function (IRF) P7SOURCE_V6. 
Event and spacecraft all-sky-survey data amassed between 
the 4 th of August 2008 and the 28 th of April 2012 (a 
Mission-Elapsed- Time interval of 239557417 to 357308125, 
equating to ~1362 days) were extracted for a source re- 
gion (SR) of 12°, centred on each BCG position (taken 
from optical data in the literature) and in the energy range 
300 MeV ^ E sC 300 GeV (thus discarding events at the 
lowest energies, where the high point spread function (PSF) 
results in considerable source confusion). A zenith cut of 
100° was applied to eliminate 7 rays from the Earth's limb, 
and all 'Source' class events were considered. 

For each candidate source the data were fitted to 
a source model constructed to describe (in addition to 
potential emission from the target BCG), the emission 
from nearby (i.e. within the SR) 2FGL sources, and 
the diffuse Galactic and extragalactic (isotropic) 7-ray 
backgrounds (thus recent models gal_2yearp7v6_v0.fits and 
iso_p7v6source.txt respectively!]] were utilised; their normali- 
sation free to vary) . The photon index and normalisation for 
objects within 1° of the candidate position were allowed to 
vary. Outside of this region but within the specified region 
of interest (ROI) of 8° , the normalisation was free to vary, 
but the photon index was fixed to the 2FGL value. Sources 
outside of the ROI but within the SR were included in the 
model, but their aforementioned parameters were fixed to 
the catalogue values. In the event that a candidate position 
coincided with a 2FGL source (within 0.2°), then that source 
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Figure 1. An all-sky map constructed from ~45 months of 300 MeV ^ E ^ 300 GeV Fermi data, scaled logarithmically. The BCG 
positions of the 114 candidates of the main sample arc indicated by the white crosses. 



was automatically removed from the model. This action pro- 
vides a useful check of the output of our analyses against the 
catalogue values for known Fermi sources, whilst allowing all 
BCGs in the sample to be treated in the same way. 

The maximum-likelihood spectral estimator GTLIKE 
was used to perform a binned fit, modelling the candidate 
source emission with a power law spectrum. GTMODEL was 
then used to obtain a model map of each candidate region 
given the result of the fit. To construct counts maps from the 
data gtbin was utilised. Binned exposure cubes were also 
generated (through use of the GTEXPCUBE2 tool) as part 
of the binned likelihood analysis chain. A binned method of 
likelihood fitting was adopted over an unbinned one, because 
it was found to be more robust, the fit converging irrespec- 
tive of how many 2FGL sources a given SR contained. 

The distribution of Test Statistic (TS) values (output 
by GTLIKE) across the sample could then be studied. A crit- 
ical TS of 25 (corresponding rou ghly to a detection signifi- 
cance of 5a (|Mattox et al.lll996t was decided upon, below 
which an upper limit on the source emission was automat- 
ically calculated. All candidate sources with TS s!16 were 
nevertheless investigated on an individual basis: TS maps 
were constructed and any peaks found thereon fitted with 
parabolic ellipsoids to localise the emission. Where appro- 
priate the data were then remodelled to constrain potential 
new offset sources. 

The ensemble of counts maps constructed from the sam- 
ple data, and that of the model maps output from fitting 
based on a null hypothesis (that is, the absence of a candi- 
date Fermi source at the position of the BCG), were stacked. 
These stacked maps could then be sliced to provide 1-D com- 
parisons between the summed model fit and the data. To 
derive an upper limit on the 7-ray flux of an average source 
in our sample using this stacking method, the stacked counts 
map was used to obtain a number of on-source counts {N on ) 



within a region of radius comparable with the PSF of the 
LAT instrument on average across the selected energy band. 
An estimate of the number of background photons in this 
region, TVb was obtained by similarly integrating the stacked 
model map. The number of off-source counts (AW) is then 
JVb I a, where alpha is taken as the ratio of the solid angle in 
the chosen radius to that of the field of view. Thus an esti- 
mate of the 7-ray signal present is given by the excess counts, 
N on — ct ■ N {{. A 95% confidence upper limit on th e signa l 
was computed following the method of iRolke et alj |2005), 
assuming only Poisson errors on the background estimate. 
Whilst a m ore sophisticated li kelihood method similar to 
that used bv lHuber et al.l l|2012T l may provide a modest sen- 
sitivity improvement, our approach is simple and robust. 

Weighted exposure maps were constructed for each can- 
didate source and stacked, to provide a summed exposure 
map from which the exposure averaged over solid angle in 
the on-region could be determined. This value was used to 
transform our signal upper limit to an upper limit on the 
7-ray flux. The rms variation in exposure across the sample 
is w6%. 

For the purpose of investigating potential weak, hard 
spectrum, candidate source emission; less contaminated by 
the diffuse background, and given that stacking is most 
promising in the signal-limited regime, where every addi- 
tional field adds potential source photons with minimal cost 
in background contamination, energy cuts above 3 GeV and 
30 GeV were applied to the data and the analyses repeated 
in an otherwise identical manner; upper limits calculated 
within regions of radii comparable to the LAT PSF at the re- 
spective lower energy bounds. Corresponding analyses were 
performed simultaneously for the control and radio-quiet 
samples. 
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4 RESULTS 

A summary of the source properties of the BCGs of the 
main sample; their respective TS values and/or upper limits 
output from the analysis procedure outlined in |3] is given 
in Table [JJ For the sake of clarity, the test statistic is shown 
as zero for sources with TS <1. An upper limit is included 
for sources with TS <25. Those BCGs with TS >25 are in: 

A 3392 

RXCJ1958.2-3011 
MACS J1347.5-1145 

As noted in the sources in A 3392 and 

RXC J1958.2— 3011 correspond to the known Fermi 
sources 2FGL J0627. 1-3528 and 2FGL J1958.4— 3012, 
respectively. The BCG in MACS J1347.5-U45 is not in 
the 2FGL catalogue. These sources, and those that fall shy 
of the TS cut, but for which TS is greater than «16 are 
discussed in £14.31 below. 
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Table 1: A summary of the 114 brightest cluster galaxies of the main sample, and their candidate Fermi 
source properties, including the upper limit on the 7-ray flux above 300 MeV. The radio emission is 
classified as core-dominated ('c'), extended- with-core-contribution ('e') or BL Lac-associated ('b'). 



Parent Cluster f (°) 6 f (°) z 1.4 GHz* Radio TS n 7-ray Flux Upper 













Flux (mjy) 




Limit (phcm 


2 - 
s 




RXCJ0000. 1+0816 


101.9 


-52.5 





04 


84 


c 





2.2 ■ 


10" 


10 


A85 


115.3 


-72.0 





056 


58 


c 





2.8- 


10- 


10 


Z235 


120.8 


-38.4 





083 


50 


c 


3.2 


5.7- 


10" 


10 


3C31 


126.9 


-30.3 





016 


1108 


e 


2.2 


5.8- 


10" 


10 


RXCJ0132.6-0804 


152.0 


-68.6 





15 


308 


c 


24.0 


1.6 


• 10 


-9 


MACSJ0242.5-2132 


206.6 


-64.1 





31 


1255 


c 





1.9- 


10" 


10 


A3112 


253.0 


-56.1 





075 


1915 


c 


11.2 


8.2- 


10" 


10 


A496 


209.7 


-36.5 





033 


121 


c 


5.8 


9.2- 


10" 


10 


RXCJ0439.0+0520 


191.4 


-26.2 





21 


82 


c 


6.0 


8.2- 


10" 


10 


NGC1650 


214.0 


-34.9 





036 


91 


c 





5.0- 


10" 


10 


S555 


243.6 


-26.3 





044 


319 


c 





2.1 • 


10" 


10 


A3378 


241.9 


-24.0 





14 


1337 


c 





2.9- 


10" 


10 


A3392 


243.5 


-20.0 





055 


4513 


b 


221.0 




- 




PKS0745-191 


236.5 


3.0 





1 


2373 


c 


13.7 


8.2- 


10" 


10 


A646 


172.7 


34.6 





13 


54 


c 





8.3- 


10" 


11 


4C+55.16 


162.3 


36.6 





24 


8284 


c 


1.6 


4.9- 


10" 


10 


Hydra-A 


243.0 


25.1 





054 


40850 


c 


2.7 


4.8- 


10" 


10 


MACSJ1133.2+5008 


150.7 


62.6 





39 


846 


e 





2.3- 


10" 


10 


A1348 


277.4 


47.1 





12 


158 


c 





1.4- 


10" 


10 


A1366 


132.6 


48.5 





12 


237 


b 





3.0- 


10" 


10 


NGC4696 


302.5 


21.6 





0098 




e 





3.4- 


10" 


10 


MACSJ1347.5-1145 


324.1 


48.8 





45 


48 


c 


43.4 




- 




A 1795 


33.9 


77.2 





062 


925 


e 





1.9- 


10" 


10 


RXCJ1350.3+0940 


344.3 


67.7 





13 


293 


c 





1.0- 


10" 


10 


A3581 


323.2 


32.9 





023 


646 


c 


7.4 


1.1 


■ 10 


-9 


A1885 


83.2 


66.6 





089 


49 


c 





1.2 ■ 


10" 


10 


S780 


341.0 


35.1 





24 


106 


c 


1.5 


3.1 • 


10" 


10 


RXCJ1504.1-0248 


355.1 


46.2 





22 


62 


c 


7.3 


9.6- 


10" 


10 


A2052 


9.5 


50.1 





035 


5500 


c 





2.4- 


10" 


10 


A2055 


8.9 


49.3 





1 


498 


b 


8.0 


8.4- 


10" 


10 


RXCJ1524.2-3154 


337.1 


20.7 





1 


50 


c 


2.0 


5.3- 


10" 


10 


RXCJ1558.3-1410 


356.6 


28.7 





097 


461 


c 


1.7 


6.8- 


10" 


10 


A2199 


63.0 


43.7 





03 


3681 


e 


1.5 


5.9- 


10" 


10 


NGC6338 


85.9 


35.4 





028 


57 


c 





3.5- 


10" 


10 


Z8193 


67.6 


34.7 





18 


134 


c 





2.7- 


10" 


10 


A2270 


83.1 


33.8 





24 


144 


c 


2.6 


5.1 • 


10" 


10 


Z8276 


58.0 


27.6 





076 


92 


c 





5.6- 


10" 


10 


RXCJ1750.2+3504 


60.6 


27.0 





17 


69 


c 





1.5- 


10" 


10 


E1821+643 


94.1 


27.4 





3 


94 


c 


1.3 


2.7- 


10" 


10 


RXCJ1832.5+6848 


99.2 


26.8 





2 


150 


c 





1.1 ■ 


10" 


10 


A3639 


347.0 


-26.3 





15 


139 


c 





2.1 • 


10" 


10 


RXCJ1931.6-3354 


5.2 


-22.5 





097 


886 


c 





2.8- 


10" 


10 


MACSJ1931.6-2634 


12.6 


-20.1 





35 


223 


c 


7.4 


1.0 


■ 10 


-9 


RXCJ1958.2-3011 


11.0 


-26.8 





12 


128 


b 


46.3 








Cyg-A 


76.3 


5.8 





056 


1590000 


e 


4.5 


7.9- 


10" 


10 


S851 


351.0 


-32.6 





0097 


139 


c 





9.3- 


10" 


10 


RXCJ2014.8-2430 


18.4 


-28.5 





16 


230 


e 





3.4- 


10" 


10 


A2390 


74.0 


-27.8 





23 


236 


c 


11.2 


4.3- 


10" 


10 


RXCJ2213.1-2754 


22.3 


-55.0 





061 


76 


c 





2.8- 


10" 


10 


A3880 


18.1 


-58.5 





058 


228 


c 


4.0 


5.7- 


10" 


10 


A2597 


65.4 


-64.9 





085 


1875 


c 





1.9- 


10" 


10 


A2627 


101.8 


-35.9 





13 


434 


b 


9.8 


8.3- 


10" 


10 


A2634 


103.6 


-33.1 





03 


1037 


e 





1.9- 


10" 


10 


RXCJ0137.2-0911 


156.3 


-69.1 





041 


178 


e 





3.1 • 


10" 


10 
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Table 1 - 


continued 










Parent Cluster 


r(°) 


&t(°) 


z 


1.4 GHz* Radio 


TS tt 


7-ray Flux Upper 










Flux (mjy) 




Limit (phcm 


V 




A262 


136.6 


-25.1 


0.017 


67 e 





A ^ . 


1 n~ 


10 


A C\ t~\ O A 

A2984 


256.6 


-68.9 


0.1 


232 c 









10 


A3017 


256.6 


-65.7 


0.22 


127 c 


3.3 


1 • 

<J. 1 


1 fl~ 


10 


Z808 


175.4 


-47.3 


0.17 


393 e 





Q fi • 


1 fl~ 

J-U 


11 


A407 


150.7 


-19.9 


0.048 


661 c 





1 7 • 
1 . 1 


1 fl~ 

J-U 


10 


RXCJ0331. 1-2100 


212.3 


-53.2 


0.19 


168 c 





9 t\ ■ 


1 n~ 

J-U 


10 


NGC1399 


236.8 


-53.6 


0.0051 


2500 e 


2.9 


A 8 . 


JU 


10 


A3165 


226.4 


-51.4 


0.14 


921 c 





1.4 • 


1 n~ 


10 


RXCJ0359.1-0319 


193.4 


-39.3 


0.12 


188 e 





9 8 • 


1 fl~ 

JU 


10 


RXCJ0425. 8-0833 


203.4 


-36.2 


0.04 


112 e 





9 ^ ■ 


1 fl~ 

JU 


10 


MACS J0429. 6-0253 


197.9 


-32.6 


0.4 


139 c 





1.0 


1 n~ 

JU 


10 


A499 


218.5 


-38.3 


0.15 


158 e 





A 8 . 


1 n~ 

JU 


10 


RXCJ0501.4+0H0 


198.5 


-23.7 


0.12 


121 c 





^ 9 • 


1 n~ 

JU 


10 


MACSJ0520.7-1328 


215.3 


-26.1 


0.34 


93 c 





9 8 . 


JU 


10 


A3360** 


249.4 


-30.7 


0.085 


222 c 





2.2 • 


1 fl~ 

JU 


10 


S547 


254.2 


-30.0 


0.051 


143 e 





1 • 


1 fl~ 

JU 


10 


A3380 


257.2 


-27.3 


0.055 


1455 e 


1.0 


^ 1 • 


JU 


10 


A3396 


250.0 


-21.6 


0.18 


614 e 


2.5 


2.4 • 


1 n~ 

JU 


10 


A795 


217.1 


40.2 


0.14 


116 c 


1.1 


1.0 


JU 


10 


S617 


253.3 


23.3 


0.034 


105 e 


7.3 


7 • 


1 fl~ 

JU 


10 


A907 


249.4 


33.3 


0.17 


69 c 







1 fl~ 

JU 


10 


A923 


205.3 


53.3 


0.12 


73 c 





1 Q • 


1 fl~ 

JU 


10 


Z3179 


228.7 


53.1 


0.14 


94 c 


6.6 


7 A ■ 


1 n~ 

JU 


10 


Z3916 


144.7 


54.4 


0.21 


66 c 





q n . 


JU 


11 


A1361 


153.4 


66.6 


0.12 


864 e 





O.vt 


1 n~ 

JVJ 


10 


A3490 


287.8 


26.5 


0.07 


368 e 





9 n . 


1 fl~ 

JV^ 


10 


A1412 


128.4 


43.1 


0.083 


73 e 





9 n . 


1 fl~ 

JU 


10 


MACS J1206. 2-0847 


284.5 


52.4 


0.44 


161 e 





1.1 • 


1 fl~ 

JU 


10 


A1644 


305.0 


45.4 


0.047 


98 c 





iJ.O 


1 n~ 

JU 


10 


A1677 


84.0 


85.1 


0.18 


76 c 





1 7 ■ 


1 n~ 

JVJ 


10 


NGC5098 


78.8 


81.3 


0.036 


83 c 





0. 1 


JU 


10 


A3565 


313.6 


28.0 


0.012 


3000 e 





^ A ■ 


1 fl~ 

JU 


10 


TV T A /~1 n T1 ^11 o i r Ai n 

MACSJ14H.3+5212 


97.6 


60.8 


0.46 


22720 e 





3.9 • 


10~ 


10 


A2204 


21.2 


33.2 


0.15 


70 c 







1 fl~ 

JU 


10 


HerA 


23.1 


28.9 


0.15 


34000 e 


3.6 


0.0 


1 n~ 

JVJ 


10 


AT f O f n 

NGC6269 


49.1 


35.9 


0.035 


51 c 


3.0 


9 ^ ■ 


JU 


10 


RXCJ1715. 1+0309 


24.5 


22.8 


0.13 


165 c 





9 7 ■ 


1 n~ 

JVJ 


10 


RXCJ1720. 1+2637 


49.3 


30.9 


0.16 


89 e 





1.4 ■ 


1 fl~ 

JU 


10 


RXCJ1840. 6-7709 


317.3 


-25.8 


0.019 


1152 e 





9 ^ • 


1 fl~ 

JU 


10 


RXCJ1852. 1+5711 


87.1 


22.4 


0.11 


51 c 





A *\ ■ 


1 fl~ 

JU 


10 


A2312 


99.0 


24.8 


0.095 


79 c 


5.4 


6.3 • 


10~ 


10 


A3638 


355.4 


-24.0 


0.077 


258 c 





2.4 ■ 


1 n~ 

JVJ 


10 


IC4991 


359.9 


-33.2 


0.019 


109 c 


16.0 




JU 


-9 


S868 


22.9 


-29.1 


0.056 


229 c 





2.1 • 


1 fl~ 

JU 


10 


RXCJ2034. 9-2143 


23.2 


-32.0 


0.19 


82 c 





"i A ■ 


1 fl~ 

JU 


10 


RXCJ2043. 2-2144 


23.9 


-33.8 


0.2 


295 e 





1 8 . 


1 fl~ 

JU 


10 


A2331 


40.9 


-31.7 


0.079 


106 e 


5.3 


0.0 


1 n~ 

JVJ 


10 


RXCJ2101. 5-1317 


35.4 


-34.8 


0.028 


1200 e 





1.9 ■ 


10" 


10 


RBS1712 


22.2 


-38.8 


0.19 


67 c 





1.9- 


10" 


10 


APM699 


338.1 


-45.7 


0.082 


152 e 


11.7 


9.8- 


10" 


10 


RYPT91 ^1 *3 ^91 
XWVvy J Z _1_ _1_ . o-OOZ 1 


9 


AT 1 


U.Uoo 


1 uo e 




3.0- 


10" 


10 


A2384 


33.4 


-48.4 


0.096 


54 c 





3.1 • 


10" 


10 


A2389 


53.7 


-41.8 


0.15 


78 c 





3.6- 


10" 


10 


A2415 


54.0 


-45.1 


0.058 


187 c 





1.4- 


10" 


10 


A3856 


2.8 


-56.2 


0.14 


186 e 





1.6- 


10" 


10 


A2445 


83.0 


-30.8 


0.17 


77 c 


6.0 


4.3- 


10" 


10 


A3911 


336.7 


-55.4 


0.097 


418e 





1.5- 


10" 


10 
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Table 1 - 


continued 






Parent Cluster 


m 


&t(°) 


z 


1.4 GHz* Radio 


TS tt 


7-ray Flux Upper 










Flux (mjy) 




Limit (phcm -2 s _1 ) 


S1101 


348.4 


-64.8 


0.056 


472 c 





3.0 • 10" 10 


A4023 


304.7 


-31.7 


0.19 


143 e 





2.9 ■ 10" 10 


A2665 


97.0 


-53.6 


0.057 


56 c 





1.7- nn 10 



Galactic longitude 
1 Galactic latitude 

t Taken from the NVSS. If italicised, the values instead denote the 843 MHz radio flux taken from the SUMSS 
1"! Test Statistic - for clarity TS values <1 are shown as zero 

Despite a low TS, this source was not included in the stacked analysis (see i]4,2jl due to the proximity of a bright 2FGL source, as 
detailed in the text 
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MR 



10 12 



14 



in the colour scale, overlayed with contours defining the 
corresponding stacked model map for the null hypothesis. 
The maps to be summed were each centred on the can- 
didate source position (0,0), which is marked on the fi- 
nal stacked plot. The right panels show one-dimensional 
slices through both the summed counts and model maps 
(the area covered by these slices is indicated by the boxes 
drawn on the stacked maps; their width motivated in each 
case by the LAT PSF as described in f£3 to match the re- 
gions within which upper limits are computed: 1°, 0.4° and 
0.2° for the ^300 MeV, ^3 GeV and ^30 GeV datasets, 
respectively). Sources with a high individual TS (>16), and 
those within 0.2° of a catalogue source are not included in 
the stack. The field containing the BCG A 3360 has also 
been removed, due to the proximity (within 1°) of the 
very bright ((3.71 ± 0.07) x 10~ 8 phcm~ 2 s~ 1 ) 2FGL source 
2FGL J0538.8— 4405. The model provides a reasonable de- 
scription of the data, and there is no evidence for an excess 
close to the origin. Applying the method described in Sj3[ 
the derived 95% upper limits on the 7-ray flux of an av- 
erage source within the class of BCGs, given the cuts ap- 
plied and within a region of size comparable with the LAT 



PSF, are found to be 7.4x 10 _11 phcm" 



(within a 1° ra- 



dius for 300 MeV < E < 300GeV), 7.1x lO-^phcm^s -1 
(within a 0.4° radius for 3 GeV < E sC 300 GeV) and 



TS 



S^xlO-^phcm" 
E < 300 GeV). 



(within a 0.2° radius for 30 GeV < 



Figure 2. Histograms showing the comparison between the dis- 
tribution of detection significance across the samples described 
in [T2] The main sample of BCGs (black line) is included on both 
plots. Top Panel: The (averaged) control sample of random posi- 
tions on the sky (dashed violet line). Bottom Panel: The scaled 
extended cluster sample (dot-dashed green line) and radio-quiet 
cluster sample (dashed blue line). The (two) sample BCGs with 
existing Fermi detections have been added in red. 

4.1 Distribution of Detection Significance 

The distribution of \/TS (corresponding approximately to 
the statistical significance, see Sj3| for each sample analysed 
is shown in Figure [2] The extended and radio-quiet clus- 
ter samples have been scaled appropriately, and the control 
samples averaged, so that comparisons may easily be made 
with the main sample of 114 objects. Common to all distri- 
butions is a clear peak at 0, and also a tail containing several 
signals at up to ~3cr. However this tail is more pronounced 
for the main sample - implying that a stacking analysis may 
be worthwhile (see H4.2I below) - and outliers are present: 
the high-TS sources are visible above 4<r (see i)4.3l below). in- 
cluding the two Fermi sources shown in red; all members of 
the main sample. The only exception is the extended sample 
BCG in S 753, for which a TS of 19.9 was derived. Inspection 
of the field containing this object reveals a highly confused 
region of diffuse emission, likely associated with the Galactic 
background, and as such this result is not believed. 

4.2 Stacking Analysis 

The stacked analyses of the main sample above the three 
energy cuts - described in Sj3] - are illustrated in Figure [3] 
In each case the left panel shows the stacked counts map 



4.3 Sources of Note 

4.3.1 2FGL Sources 

The Fermi source 2FGL J0627. 1—3528 is associated with 
the BCG in A 3392, which is classified as a BL Lac ob- 
ject. It is listed with a 24-month detection significance of 
12.6, 100 MeV < E < 100 GeV flux of (1.50 ± 0.18) x 
10" 9 phcm- 2 s _1 and photon index of 1.93 ± 0.09. With 21 
months more data (and in an energy band 300 MeV ^ E ^ 
300 GeV) we calculate an increased detection significance 
of 14.9, and a consistent (within errors) photon index, as 
expected. 

The Fermi source 2FGL J1958.4— 3012 is associated 
with the BL Lac object RXC J1958. 2-3011, and may be 
a BCG (see § EU- It is listed with a 24-month detection 
significance of 5.1, 100 MeV < E sC 100 GeV flux of 
(5.94 ± 1.78) x 10 _10 phcm~ 2 s~ 1 and photon index of 1.9 
± 0.2. With 21 months more data (and in an energy band 
300 MeV ^ E ^ 300 GeV) we calculate an increased de- 
tection significance of 6.8, and a consistent (within errors) 
photon index, as expected. 

4-3.2 Offset Detections: Possible Counterparts 

TS maps were constructed for candidate sources with 
TS ~16, to better localise the source position. These can 
be seen in Figure H for the BCGs in MACS J1347.5-1145 
(TS = 43.4), RXCJ0132.6-0804 (TS = 24.0) and IC4991 
(TS = 16.0), and illustrate the effect of neglecting to add 
the candidate to the model: significant sources are revealed 
and localised where the baseline model (containing all local 
2FGL sources, and the diffuse Galactic and extragalactic 
backgrounds, see $3} fails to describe the 7-ray emission. 
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Figure 3. Stacked counts maps for the main sample of BCGs above an energy (top) 300 MeV, (middle) 3 GeV and (bottom) 30 GeV, 
overlayed with contours representing the corresponding stacked model maps (for the null hypothesis). The white cross indicates the BCG 
position in all fields. The stacked maps have been smoothed with a Gaussian kernel of rms 0.2°. A cut in candidate source TS above 16 
has been applied, and additionally a proximity cut of 0.2° (to remove fields where the candidate position is coincident with a Fermi-hAT 
source). Neglecting also fields for which a 2FGL source brighter than 1 X 10 -8 phcm - 2 s — 1 is found to be within 1° of the central position 
(see !|4.2t . these stacks then comprise 109 (top) and 111 (middle; bottom) individual fields. The boxed regions indicate the area of 
the slices (to the right of each stack) through both maps in the horizontal and vertical directions, allowing a 1-D comparison between 
the data (black points) and null model (blue line). 
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These peaks in TS were fitted with parabolic ellipsoids and 
the best-fit peak positions are illustrated in the figure. The 
contours represent the 95% confidence regions for the posi- 
tions of the sources. 

Adding a candidate at the position of the BCG in 
MACS J1347.5-1145 (in the energy band 300 MeV < E < 
300 GeV) yields a photon index of 2.76 ± 0.01 and a flux 
of (3.81 ± 0.70) x 10" 9 phcm" 2 s"\ The significant source 
in the field however, is clearly offset from the BCG posi- 
tion (the fitted peak in TS is located at Right Ascension 
(RA) = 207.33°, Declination (Dec.) = -11.55°) and so un- 
likely to be associated. Consulting the literature reveals the 
Fermi-LAT detection of a GeV flare consistent (within the 
quote d statistical errors) with this location in November 
2011 (jATel #37881 ). and the flat-spec trum z = 0.34 Quas i 
Stellar Object (QSO) PKS 1346-112 (| Jackson et al.ll2002h . 
is given as a possible counterpart. Following this detection, a 
possible Wide-field Infrared Survey Explorer (WISE) blazar 
counterpart w as identified, coincident with the radio source 
jATel #40861 ). Re -running the likelihood analysis with this 
new source included in the model yields a negligible TS at 
the BCG position, and a TS of 179.6 (corresponding to a 
detection significance of «13.4<r) at the peak position, with 
a corresponding photon index of 2.44 ± 0.03 and flux of 
(3.04 ±0.15) x 10^ 8 phcm~ 2 s~ 1 . Despite the high TS value, 
the transient nature of the source may account for its ab- 
sence from the 2FGL. 

Adding a candidate at the position of the BCG in 
RXC J0132. 6-0804 (in the energy band 300 MeV < E < 
300 GeV) yields for the source a photon index of 2.20 ± 0.02 
and a flux of (1.42 ± 0.50) x 10~ 9 ph cm s . There appear 
to be two significant sources in the field, the more central of 
which (the fitted TS peak is located at RA = 23.25°, Dec. 
= —7.98°), is 8' from the BCG and may be associated with 
a member of the parent cluster. The most significant peak 
is fitted to a position at RA = 21.78°, Dec. = -8.37°. The 
BL Lac FBQS J012 7-0821 at a redshif t of 0.36 is a plausible 
7-ray counterpart (IPlotkin et al.l l2008), A re-analysis of the 
field using a model containing these two localised sources 
nullifies the BCG TS, and results in a TS of 28.0 for the 
closer peak, and 65.0 for the remaining one, corresponding 
to «5.3cr and «8cr detections, respectively. A photon index 
of 2.36 ± 0.03 and flux of (1.14±0.06) x 10" 8 phcm'V 1 is 
calculated for the latter, and for the former a photon index 
of 1.98 ± 0.03, and a flux of (3.06±0.24) x 10~ 9 phcm~ 2 s~ 1 . 

Adding a candidate at the position of the BCG in 
IC4991 (in the energy band 300 MeV < E < 300 GeV) 
yields a photon index of 2.379 ± 0.009 and a flux of 
(1.64±0.62) x 10" 9 ph cm s . There appear to be two po- 
tential non-2FGL sources of 7 radiation in the field. For the 
first, the BCG lies just outside the 95% confidence contour 
(the fitted peak in TS is located at RA = 304.30°, Dec. = 
—41.22°). If this TS peak represents a genuine 7-ray source, 
it may be associated with the candidate, or perhaps with 
the ROSAT source 1RXS 201731.2-4 11452, a BL Lac object 
6' removed l|Kollatschnv et aljliooj ). The second TS peak 
(fit position: RA = 307.27°, Dec. = -42.72°) may indicate 
7-ray emission associated with an 81 mjy S UMSS source 
5' rem oved, that is not clearly identified (see iMauch et al.l 
(2003)). Remodelling the data to account for these two 
sources renders the candidate BCG emission insignificant 
(reducing the TS to 2.7), and results in a TS of 13.9 for the 



Table 2. A summary of non-2FGL Fermi detections of sources 
in fields within our sample, giving the fitted position (and error) 
of the TS peak, and the output TS after remodelling the data. 
No attempt is made to account for trials factors due to the large 
sample considered, but we note that a TS value of >25 still cor- 
responds to a post-trials significance of >4o\ 



Name TS* RA 1 " (°) Dec* (°) Error (°) 



FERMI J1350-1140 1 " 1 " 


179.8 


207.33 


-11.55 


0.11 


FGLJ0133.0-0758 


28.2 


23.25 


-7.98 


0.16 


FGL J0127. 1-0822 


64.8 


21.78 


-8.37 


0.14 


FGL J2029.0-4243 


41.3 


307.27 


-42.72 


0.22 



* Test Statistic 

* Right Ascension 

* Declination 

^ Detected by the Fermi Large Area Telescope Collaboration in November 
2011 (ATcl #3788) 



closer peak, and 41.1 for the peak revealed toward the edge 
of the field, roughly corresponding to statistical significances 
of 3.7cr and 6.4cr, respectively. A photon index of 2.2 ± 0.2 
and flux of (8.58 ±3.15) x 10" 9 ph cm 2 s 1 is calculated for 
the significant source. 

Table [2] summarises the detections confirmed after a 
remodelling to account for sources of 7-ray emission in the 
three fields described above was carried out. The TS maps 
were re-calculated based on the new output models, which 
were found to describe the emission adequately. 

4-3.3 j-ray Emission from A 2055? 

Using the 3 GeV < E < 300 GeV data, analysis of the BCG 
in A 2055 (which is classified as a BL Lac object) results in a 
TS of 15.2. The corresponding TS map is shown in Figure[5] 
A hint of emission from the position of the BCG can be 
clearly seen. The position of the BCG is contained within 
the 95% confidence region for the source location. With this 
candidate alone added, the model sufficiently describes the 
7-ray emission within the field. The source is not detected in 
the energy band 30 GeV < E s£ 300 GeV, and the TS drops 
to 8.0 using the larger 300 MeV < E sC 300 GeV dataset. 



5 DISCUSSION 

There is no statistically significant evidence for emission 
within the stacked sample. The upper limits derived on an 
average source within the class of BCGs represent at least 
an order of magnitude improvement on the individual upper 
limits presented in Table [T] The derived limits are consis- 
tent with the expectation of a ~ 1/ WN^bj scaling for the 
background-limited lower energy bands and ~ l/iV bj for 
highest energies where the Fermi-LAT becomes signal lim- 
ited even in the stacked sample, where iV D bj is the number 
of fields stacked. 

To derive a limit on the 7-ray to radio flux ratio 
(FiGev/Fi.4GHz) for the core emission of the sample we con- 
sider the core radio flux (so me fraction of the e mission given 
in Table [1} of each BCG (jHogan et alJlprepn . and find an 
average radio flux of 2.37 x 10 -15 ergcm~ 2 s _1 and an aver- 
age luminosity of 1.8 x 10 41 ergs _1 (at an average redshift of 
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Figure 4. TS map for candidate sources in MACS J1347.5— 1145 
(top), RXCJ0132. 6-0804 (middle), and IC4991 (bottom) pro- 
duced from 300 MeV ^E^ 300 GeV data. The BCG position is 
indicated by the central black/grey cross. Peaks in TS have been 
fitted with parabolic ellipsoids, the fit position for which is shown 
by a white triangle and given in table [2] The 95% confidence 
contours on source locations are shown in white. 



8.5 
8 
7.5 

7 

§ 65 

-E 6 
o 

CD 

° 5.5 
5 
4.5 
4 

232 231 230 229 228 

Right Ascension (°) 

Figure 5. TS map for the candidate source in A 2055, produced 
from 3 GeV s: E ^ 300 GeV data. The BCG position is indicated 
by the central black/grey cross. The 95% confidence contour on 
the 7-ray source position is shown in white. 

0.13). Given the upper limit on the stacked emission, the av- 
erage 7-ray flux evaluated at 1 GeV and assuming a photon 
index of 2, is less than 3.55 x 10 -14 erg cm -2 s _1 , and the 
average 7-ray luminosity is less than 1.35 x 10 42 ergs _1 . The 
average -F\gcv/-Fi.4GHz of the sample is therefore <15. The 
7-ray flux-limit for the sample corresponds to only 0.1% of 
the flux of NGC 1275, which has FioeV /Fiagkz -120 (~ 90 
excluding 7-ray flares). The average core radio flux is «1% 
of that of NGC 1275. The lack of a 7-ray signal in the stack 
therefore implies that the average core-radio selected BCGs 
have spectral energy distributions with a different shape to 
that of NGC 1275, with significantly lower 7-ray flux relative 
to core radio flux, though the current highly variable state 
of 7-ray emission in this active galaxy (|Kataoka et aljfeoiol ) 
may account for the difference in Figcv/^i^ghz. 

NGC 1275 and M 87 are the closest cooling-core BCGs, 
so the detected 7-ray emission is unlikely to be beamed 
toward us if they are drawn from a random distribution 
of sight lines. However, a role for beaming in NGC 1275 
cannot be excluded and may play a role in its unusual 
7-ray brightness (see above). It has been postulated (see 
iFalceto-Goncalves et ail (|2010T )) that the central AGN com- 
prises a viscous accretion disk tilted with respect to the ro- 
tating super-massive black hole, bringing about precessing 
jets as a result of torques acting on the disk. Very Long 
Baseline Interferometry (VLBI) maps of the radio source 
can place constraints on the inclination of the system: to 
be satisfied, some curvature of the .j et away from the line 
of sight is implied dDunn et al.l 1200(1 ) . This said, for those 
BCGs with bona fide BL Lac nuclei, for example A 2055 and 
A 3392, which are significantly more distant and comparable 
to most of the parent sample, the probability of observing a 
few sources close to the jet direction is significant. 

We can assess the number of local clusters hosting a 
Fermi-LAT AGN by cross-correlating the Abell cluster cat- 
alogue of 5250 clusters (including the supplementary sys- 
tems) l| Abell et all 1 19891 ), with the 2FGL. Of the 85 coin- 
cidences within 20', only 11 are consistent with the 2FGL 
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source being a cluster member and of those only four be- 
ing the BCG (NGC1275 in A426/Perseus, PKS 0625-35 
in A 3392, PKS 2035-714 in A 3701 and PKS 2316-423 in 
Sllll). The latter two clusters are not detected above the 
X-ray flux limit for the eBCS or REFLEX catalogues, so 
were not considered in the analysis above. Therefore, the 
chance of finding a Fermi source in any particular cluster is 
relatively small but is significant. 

The added complication with the analysis is that the 
Fermi-LAT spatial resolution and source density (~0.05 per 
square degree) implies chance coincidences on the order of 
a few when more than 100 positions are considered. This 
is illustrated in Figure 3J given that the TS output from 
the initial binned likelihood analysis of the field was in each 
case not exclusively associated with the BCG under scrutiny, 
but a nearby source of emission not described by the model. 
The angular resolution of the Fermi-LAT is sufficient (par- 
ticularly above 1 GeV) to isolate nearby background sources 
from the position of interest in most cases, but the likelihood 
analysis is not robust enough to accurately disentangle the 
7 emissions in such an event. 

Four detections (that is, of TS>25 emission) have been 
made of sources not contained in the 2FGL (see Table [2} . 
For those above a TS of ~50, their absence from the 2FGL 
implies temporal variability in the emission, as appears to be 
the case for the source FERMI J1350— 1140. A more detailed 
study of these sources is left to future work. The tentative 
evidence for emission from the BCG in A 2055 will be refuted 
or validated in time: if the significance of emission from the 
source position continues to increase, it will be the fourth 
such source to be detected using Fermi; the second BCG 
hosting a BL Lac. 



6 CONCLUDING REMARKS 

The Fermi-LAT detection of a few sources in clusters of 
galaxies sets important limits on the energetics and emis- 
sion mechanisms at work in the cores of the most massive 
member galaxies. The detections of BCGs using Fermi are 
strongly biased toward the most radio bright systems or 
those that are dominated by a beamed core. Our statistical 
upper limit from a larger sample of sources with lower radio 
flux indicates that the average core-radio selected BCG is 
more than an order of magnitude less 7-ray luminous than, 
for example, NGC 1275. Furthermore the ratio of 7-ray to 
radio flux is <15, compared with -FiGev/^i^GHz ~120 for 
NGC 1275. The BCG in A 2055, which hosts a BL Lac, may 
soon be detectable in HE 7 rays (given a longer exposure 
using Fermi ), and if so, will corroborate this conclusion. 

Detections of a number of new sources are declared; 
though none of these appear to be unambiguously associated 
with the BCG, the slightly-offset emission in the field of 
RXC J0132.6— 0804 might be associated with the candidate, 
or else another cluster member galaxy, and in the field of 
IC4991 may (if genuine) in part originate from the BCG. 
This work also illustrates then, the need to carefully assess 
all the possible sources of emission that could explain the 
Fermi-LAT detections found in clusters before attributing 
them to (for example) cluster merger-related shocks or DM 
annihilation. 
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